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We develop a new projected wave function approach which is based on projection operators in 
the form of matrix-product operators (MPOs). Our approach allows to variationally improve the 
short range entanglement of a given trial wave function by optimizing the matrix elements of the 
MPOs while the long range entanglement is contained in the initial guess of the wave function. The 
optimization is performed using standard variational Monte Carlo techniques. We demonstrate the 
efficiency of our approach by considering a one-dimension model of interacting spinless fermions. 
In addition, we indicate how to generalize this approach to higher dimensions using projection 
operators which are based on tensor products. 
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One of the basic problems in physics is to find the 
ground state of strongly correlated electron systems. 
These are of central importance in several areas of science 
and technology, including solid state physics, quantum 
chemistry as well as nanotechnology. Analytical exact so- 
lutions are only available in few cases and one thus often 
relies on numerical studies. Efficient numerical methods 
also exist only for special cases. Quantum Monte Carlo 
(QMC) simulations [l[ are for instance only efficient in 
the absence of any negative sign problem which occurs 
in most fermionic and frustrated systems. By contrast, 
exact diagonalizations are applicable in general fermionic 
models. However, the size of the corresponding Hilbert 
space grows exponentially with the system size and con- 
sequently only very small systems can be considered. In 
one dimensional systems, very accurate results have been 
obtained using the density matrix renormalization group 
(DMRG) @. The DMRG algorithm is closely related to 
the concept of matrix-product states (MPS) [!, 0] which 
can be generalized to higher dimensions using tensor- 
product states (TPS) @. Different TPS based or related 
approaches have been introduced recently, e.g., the pro- 
jected entangled pair states (PEPS) d, 0], string-bond 
states 0, scale- renormalized MPS j|, tensor entangle- 
ment renormalizatio n g roup (TRG) (lol H]| , correlator- 
product states [12], LL3( and entangled-plaquette states 
[14H16I ] . It has been shown that Monte Carlo sampling in 
combination with MPS or TPS yields an efficient method 
to simulate spin systems 17|, [l8[. Although there ex- 



of matrix-product operators (MPOs) [2ll - l23j which are 
the operator analogue of MPS. The quality of the trial 
wave function can then be successively improved by in- 
creasing the dimension of the matrices used in the MPO. 
We will refer to these states as matrix-product projected 
states (MPPS). Expectation values of physical observ- 
ables are computed using Monte Carlo sampling. We 
use a variational Monte Carlo (VMC) scheme to opti- 
mize the MPO by minimizing the energy. The proposed 
MPPS approach can be generalized to higher dimensions 
by using tensor-product based projection-operators. Fur- 
thermore, the fermionic structure of the wave function 
can be incorporated by using, e.g., a simple Slater deter- 
minant constructed from mean-field Hamiltonians. This 
approach can work for frustrated or fermionic systems 
which cannot be dealt with using QMC. We discuss the 
approach for simplicity for the case of a one-dimensional 
system in detail and then indicate how to generalize it to 
higher dimensions. 

The MPS representation of a quantum state on a chain 
of length L with periodic boundary condition is given by 
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where A!£, are \ x X matrices. The physical index jy. = 
1 . . . d represents d local states at site k. This state can 
be formally rewritten as 



ist applications to fermionic systems as well, the large 
amount of entanglement in fermionic wave functions is 
very difficult to capture, see e.g., Refs. 1|1 20j |. 

In this Letter, we propose a novel construction of trial 
wave functions to efficiently represent ground states of 
correlated spin and electronic systems. Our starting 
point is a known wave function which we use as a first ap- 
proximation of the ground state, e.g., a mean- field wave 
function or any other better guess we might have. We 
then introduce projection operators represented in terms 
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in terms of the projector MPOs, P?A (= A 3 ^ \jk)(jk\), and 



6g roduct ) is a site-factorized (product) state defined as 
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FIG. 1: (Color online) Schematic representation of a) a 
matrix-product state (MPS) and b) a matrix-product pro- 
jected state (MPPS) with periodic boundary condition in a 
one-dimensional system. |0o) is a guess to approximate for 
the trial wave function \ip). 



Clearly, any site-factorized state can be represented us- 
ing x = 1 matrices which simply renormalize the weight 
for local states on each site individually. The entire en- 
tanglement in |?/; MPS ) is expressed by the matrices Al^ 
and thus a more entangled state requires a larger x- The 
bond dimension x needed scales exponentially with the 
entanglement entropy S, with S(= Trp rc( j log p r ed) being 
the von-Neumann entropy of the reduced density matrix 
Pred- The success of the MPS representation in one di- 
mensional systems is based on the fact that the entangle- 
ment in ground-state wave functions is low and thus they 
can usually be well approximated using a small bond di- 
mension x 01 • The main idea of this letter is to choose 
an initial wave function |</> ) which is already a better ap- 
proximation of the ground state than the site-factorized 
state |0Q roduct ) in Eq. @. The improvement over an 
MPS representation is that a part of the entanglement is 
contained in \(/)o}. In other words, the projector MPOs 
only require modifying the local entanglement since the 
long range entanglement is already included in \<fio). This 
yields a much faster convergence in terms of the bond 
dimension compared to MPS (see the example discussed 
below). Furthermore, it allows us to express a larger class 
of states than MPS. In particular states which violate the 
area law, e.g., critical states, can be represented by using 
the MPPS approach. A pictorial representation which 
compares the MPS with the MPPS is given in Fig. [TJ 

As a specific example, we now discuss a model of 
spinless fermions at half-filling and show that the trial 
wave function is substantially improved over the stan- 
dard MPS approach by choosing |</> ) to be a mean-field 
Slater determinant. The Hamiltonian is given by 



H = -t 
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where hi — c\c i and V, t > 0. Using the Jordan- 
Wigner transformation, the model is equivalent to 
the XXZ spin- 1/2 Hamiltonian on the chain H = 
£(«) [tf + SiS+ +1 ) + J -i{sl + 1) • (sf +1 + 1)] 



where s- = 
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arc the Pauli matrices and arc 



related to the fermions by sf = 2rii — 1, sf — cj, and 
= Cj. Thus J z = V and J± = —2t. The discussion 
below will be largely in terms of the fermionic model, 
although all the results are valid for the spin system as 
well. 

In the non-interacting case (i.e., V = 0), the Hamilto- 
nian (j4]) has a simple ground-state wave function which 
we choose as the initial wave function for the MPPS: 
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Here tp denotes the Fermi energy, is the dispersion 
relation and the local states are given by = 0,1. 
Sj 1 j L is the Slater determinant for the configuration 
. . . Using this initial wave function, the MPPS 
Ansatz reads 
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There are two points in parameter space at which the 
ground-state wave function can be represented exactly 
with x = 1- Firstly, if V — 0, the ground state is given 
by |^o latcr } an d thus does not require any modifications 
(i.e., j4j£] = 1 for all k and jk)- Secondly, if V — > 00, 
the ground state is a two-fold degenerate charge density 
wave (CDW) state (which is a product state in this limit). 
This state is simply obtained by choosing for example 
A[2k] = ^pfc+i] = 1 f° r au ^ ( au other elements being 
zero) - everything except one of the two CDW states is 
projected out from the Slater determinant. Note that in 
the former case, the state has algebraic correlations and 
thus an MPS requires a bond dimension which roughly 
scales linearly with the system size [241 ] . In the latter 
case, the MPS can also represent the ground state with 
X=l- 

Additionally, we consider for comparison another pop- 
ular variational state with density-density Jastrow cor- 
relation factor. The Jastrow-type wave function can be 
defined as 



Jastrow 



} = cxp 
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with r\ij — In frg 



Here r,-j = |sin (j^(xi — Xj))| 

and n' 1 = 1 — hi. The two parameters vp with the nearest 
and second-nearest neighbor indices (3 are for short range 
hole- hole repulsion if these values are less than 1. The 
factor r°j is for long range correlations and repulsive if a 
is positive. In fact, we could use the optimized |^i Jastrow ) 



3 



0.14 
0.12 
=0.10 
0.08 
0.06 
0.04 
0.02 
0.00 

10 1 



W 10 



-■ — Jastrow 
-•-MPS(z=2) 

-o-MPS ( Z =3) 
-•■■-MPPS ( Z =2) 
-a- MPPS 0c=3) % 



a) 




; v/t 



6d 



10 



10 



v/t=i.o b) 

— •— MPPS — O— MPS 
V/t=2.0 

—■—MPPS — □— MPS 



-EX 

-CL CL 



° °-Oc 



— i 1 r 



8 9 10 



X 



FIG. 2: (Color online) a) Comparison of the optimized ener- 
gies for different trial states with exact diagonalization results 
for the t — V model on a L — 18 chain with periodic boundary 
condition, b) Comparison of the optimized energies with infi- 
nite time-evolving block decimation results for different bond 
dimensions \ using the MPS and the MPPS for V/t — 1, 2 on 
a L = 202 chain. Eq is the exact ground-state energy. 



act ground-state energies are obtained using sparse ma- 
trix diagonalization. The advantage of the MPPS ap- 
proach over the regular MPS is especially clear in the 
regime of small V/t. Here, the V — ground state used in 
the MPPS is a much better approximation of the ground 
state compared to the site-factorized wave function used 
in the MPS. Thus it pays off that part of the entangle- 
ment is already contained in |0Q lator ) making the MPPS 
more efficient. The MPPS yields lower energies in the 
regime of large V/t as well, even though the advantage 
is less pronounced. The simple density-density repulsive 
Jastrow-type wave function is a good approximation of 
the ground states for small V/t, however, it fails to de- 
scribe the gapped phases in the regime V/t > 2. fn par- 
ticular, it does not capture the phase transition to the 
charge ordered state at V/t = 2. Figure [5b shows the in- 
dependence of the energies and compares the results to 
the exact ground-state energies E in the thermodynamic 
limit. The energies for infinite systems are obtained us- 
ing a infinite time-evolving block decimation (iTEBD) 
[27| algorithm using MPS with large bond dimensions 
(x > 300). fn the gapless phase (V/t = 1), we find that 
the MPPS results show much faster convergence to the 
ground state than the MPS. The MPPS with \ = 5 al- 
ready has an energy which is very close to the ground 
state (error off I0~ 4 ) while the MPS displays compara- 
bly slow convergence up to \ = 10. Even at the critical 
point (V/t = 2), the MPPS gives significantly better re- 
sults than the MPS. 

In order to further test the quality of the wave function, 
we compute the density-density correlation function 
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as the initial wave function for the MPPS approach (in- 
stead of |<^o lator })- However, here we use |.(/, Jastrow } only 
as a benchmark to compare to the Slater determinant 
based MPPS. 

We now use the VMC method to optimize the en- 
ergy using the three different trial wave functions (MPS, 
MPPS and Jastrow) and to compute their correlation 
functions as well. The energy optimization is based on 
the stochastic reconfiguration (SR) method developed in 
order to optimize many parameters [25| . To stabilize the 
SR method, we also use a truncation technique for irrele- 
vant variatio nal p arameters [26j and a stochastic anneal- 
ing approach 17| . The scaling of the numerical complex- 
ity resulting from calculating the matrix product scales 
as x 3 - Note that all the matrix elements in A?u of the 
MPPS could be the variational parameters to be opti- 
mized. However, for simplicity, we assume the matrices 
to be real, symmetric and translationally invariant here. 

Figure [5^ compares the optimized energies for Hamil- 
tonian Q using a L = I 8 chain with periodic boundary 
condition with the exact ground-state energy Eq. The ex- 



Figure [3] shows the obtained correlation function C nn for 
the MPS and the MPPS using the same bond dimen- 
sion with the exact results. The exact results are cal- 
culated using the iTEBD algorithms with large bond di- 
mensions (x > 300). Note that to avoid the staggered 
decay of the correlation function arising from the 2kp 
components, we only show the correlation function with 
odd sites. Let us first discuss the gapless regime shown 
in Fig. In the V — case, the MPPS results are 
exact (up to the statistical errors) and we obtain the ex- 
pected R~ 2 decay of the correlations. For finite V/t, the 
exponents change and the correlation function is analyt- 
ically given by C nn (R) = + (-1) r ^k where A and B 
are non-universal constants and kT 1 — — arccos (— V/2t) 
[28j |. The MPPS results agree with the exact algebraic 
decay over ~50 lattice sites while the regular MPS cap- 
tures only about ~5 lattice sites (using the same bond 
dimension). In order to understand the fast convergence 
for the MPPS in Figure [2Jd, it is interesting to compare 
the behavior of the MPS and the MPPS at long distances 
in Figure^. The MPS contains no long range entangle- 
ment and thus the correlation functions saturates to a 
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FIG. 3: (Color online) Density-density correlation function 
C nn (R) for the MPPS and the MPS with the same bond di- 
mension as a function of distance R in a) gapless and b) CDW 
regions on a L = 202 chain. The green lines are obtained from 
the MPS with a large bond dimensions (x > 300) on infinite 
chains. For the chain with a finite length, the periodic bound- 
ary condition has been considered in the distance R. 



constant as R goes to infinity. In the MPPS, we start 
from a free fermion wave function in which the local en- 
tanglement is modified by the matrix-product projection- 
operators. Thus the correlation functions fall of oc R~ 2 at 
long distances. This shows us intuitively why the MPPS 
converges faster in terms of the bond dimension The 
matrices do not need to transform a product state into 
one with algebraic correlation but they only need to mod- 
ify the exponent of a wave function which is already crit- 
ical. In Fig. we also present the correlation function 
in the gapped CDW region. The CDW order parame- 
ters, enhanced as increasing V/t, can be estimated from 
the converged values of the correlation functions at large 
distance. In addition to the resemblance on the corre- 
lation functions between the MPPS and exact results, 
the MPPS also obtains a better estimation for the order 
parameter than the MPS with the same bond dimension. 

The natural extension of the MPPS to two dimensions 
is to use tensor-product projected states (TPPS). The 
trial wave function is then formally written identical to 
Eq. ([5J with the amplitudes Wj 1 j 2y ...j N of the projections 
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Here TrjS are rank z-tensors with z being the coordina- 
tion number of the lattice and N lattice size. The tenso- 
rial trace tTr is a generalization of the matrix-trace and 
implies a contraction over all virtual bond indices. The 
main difficulty is that the tensorial trace cannot be ex- 
actly evaluated anymore. This is because the numerical 
effort to contract the tensors grows exponentially with 
the linear dimension of the systems. To overcome this 
problem, we can utilize several approximate ways to cal- 
culate the contraction pSll f32| . One way to approximate 
the trace is to use the tensor entanglement renormaliza- 
tion group method 32]. The numerical complexity de- 
pends then on the lattice, e.g., applied to a honeycomb 
lattice, and that of the approximate evaluation of the 
tensorial trace scales as x 6 - The generalization of the 
method to two dimensions is now being studied 33 1. 

In conclusion, we have introduced a projected wave 
function approach which is based on matrix-product 
projection-operators. The approach is as a generaliza- 
tion of the MPS in which the initial state is not neces- 
sarily assumed to be close to a site-factorized state. We 
have demonstrated that the quality of the wave functions, 
MPPS, is superior to the MPS with the same bond di- 
mension by comparing energies and correlation functions 
for a one-dimension model of spinless fermions. The ap- 
proach can be generalized to higher dimensions by using 
tensor-product based projection-operators. Here we ex- 
pect our approach to be particularly useful because the 
computational complexity grows very quickly with the 
dimension of the tensors and thus it is very beneficial 
to express part of the entanglement by using a suitable 
initial wave function. 
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